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ABSTRACT 

We present an analysis of the Lya forest toward 3C 273 from the Space Telescope 
Imaging Spectrograph at ^ 7 km s^^ resolution, along with re-processed data from 
the Far Ultraviolet Spectroscopic Explorer. The high UV flux of 3C 273 allows us 
to probe the weak, low z absorbers. The main sample consists of 21 H I absorbers 
that we could discriminate to a sensitivity of log Nhi ~ 12.5. The redshift density 
for absorbers with 13.1 < log TV/// < 14.0 is ~ 1.5a below the mean for other lines 
of sight; for log Nhi > 12.5, it is consistent with numerical model predictions. The 
Doppler parameter distribution is consistent with other low z samples. We find no 
evidence for a break in the column density power-law distribution to log Nhi = 12.3. 
A broad Lya absorber (BLA) is within Av < 50 km s^^ and 1.3 local frame Mpc 
of two ^ 0.5L* galaxies, with an O VI absorber ^ 700 km s~^ away, similarly close 
to three galaxies and indicating overdense environments. We detect clustering on the 
Av < 1000 km s""'^ scale at 3.4(7 significance for log Nhi > 12.6, consistent with the 
level predicted from hydrodynamical simulations, and indication for a Lya forest void 
at 0.09 < z < 0.12. We find at least two components for the z — 0.0053 Virgo absorber, 
but the total Nhi column is not significantly changed. 

Key words: cosmology: observations - galaxies: intergalactic medium - quasars: 
absorption lines 



* Based on observations with (1) the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NAS 05-26555, and 
(2) the NASA-CNES/ESA Far Ultraviolet Spectroscopic Explorer 



1 INTRODUCTION 

The low redshift Lya forest offers both a challenge and an 
opportunity to study the physics of the intergalactic medium 



mission, operated by the Johns Hopkins University, supported by 
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and its role in galaxy formation and evolution. The chal- 
lenge is the requirement for space-based UV spectroscopy, 
which is undergoing a renewal with the installation of COS 
on HST. The opportunity is to observe absorption systems 
at small enough distances to be able to relate nearby galaxy 
characteristics far down the luminosity function. The cal- 
culation and interpretation of numerical mod els using the 
fluctu ating Gunn- Peterson approximation (e.g I Croft et al. I 
Il998i ) greatly benefit from constraints on the end epoch for 
absorber distribution functions such as number densities, 
column densities, Doppler parameters, clustering and rela- 
tions to galaxies. A consequence of the models is that the 
low-2 Lya forest is predicted to harbour a significant frac- 
tion of the baryons, as is a hot shock e d component (e.g . 
iDave et al. II 19991 : ICen fc Ostrikerll 19991 : lOave et al. Il200ll i. 
For a given column density, the low-z forest is predicted to 
probe larger mass overdensities than at high redshift, which 
adds complexity to comparisons of absorber property distri- 
bution functions between widely differing epochs. 

Exploration of the z < 1.6 Lya forest began in earnest 
with HST observations of the bright est QSO, 3C 273, us - 
ing the Faint Object Spectrograph (|Bahcall et al. I Il99lh 
and G oddard High Resolution Spectrograph |Morris et al. I 
Il99lh . which revealed 5 and 10 Lya absorbers, respectively. 
The main result was that there were many more low red- 
shift systems than expected from a simple extrapolation 
from z > 2 of the decline of the absorber redshift density 
over time. Concurrently, observations of the galaxy envi- 
ronments around Lya absorbers were made, again with the 
3C 273 field among the first to be probed, which revealed 
a complex relationship between absorbers and galaxies. 
The absorber-galaxy correlation was found to be stronger 
than ran dom but not a s strong as the galaxy-galaxy cor- 

'Stoc ke et al. I 
19961 : 



relation (Salzer 
1995: Lanzctta et al. 
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cent years (e.g. Chen et al. I l2005l : iProchaska et 
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Using H Lselected galaxies, however 
found the absorber-galaxy correlation even stronger than 
the galaxy autocorrelation on ~ 1 — 15 Mpc scales. 

Over the last 17 years, great efforts were made to ac- 
cumulate observations of a number of low redshift QSOs. 
A leap in observing efficiency came with the installation of 
STIS in 1997. Crucial access to higher order Lyman lines 
was provided by the launch of the Far Ultraviolet Spectro- 
scopic Explorer (FUSE) in 1999. To date, a few hundred 
low-2 Lya forest lines have been observed at resolutions 
down to ~ 7 km s~^. Such works typically involved stud- 
ies of one to a few sight lines at a time e.g. IWilliger et al. I 
(|2006l ) (Paper I), but more recently therein large analyses 
of arch i val data, including recen t work s by Lehner et al. 
('2007aVDan forth fc ShuliT (|2008l ') , and I Wakker fc Savage 
(|2009 '). The If I column density distribution shows evidence 
of a slight steepening at low z compared to that at 2: > 2 
(e.g. iLehner et al~ll2007al . and references therein), and also 
a small increase in the mean Doppler parameter (b) over 
time. An excess of broad Lya absorbers (BLAs) with b > 
40 km s~^ s-t -2^ ^ 0.5 is also observed, which is attributed to a 
larger fraction of warm-hot intergalactic medium (WHIM) 
gas at low z, consistent with model predictions; this pop- 
ulati on of BLAs is currently the subject of active study 
(e.g. iRichter. Fang fc BrvanI l2006al : iRichter et aTl l2006bl : 



iTripp et al.1l2008l : IWakker fc Savage II2009I ). Recent obser- 
vational estimates for the ba ryon contribution of the low-z 
Lya forest are 30% or more ('Penton, ShuU fc Stockel 12004 
iLehner et al. 11200781 ). [Panforth fc ShuU (2008) argue that 
low-2 O VI absorbers rev e al 10 % of the baryons in warm- 
hot gas, but iTripp et al. I (|2008l ) have shown that many of 
the O VI lines are well-aligned with narrow H I lines, and in 
those cases, the O VI absorber properties suggest an origin 
in cool, photoionise d gas. (A simila r conclusion is favored 
bv iThom fc Chen I 2008). Indeed, lOppenheimer fc Dave I 
|2009h have argued that most low-z O VI systems are pho- 
toionised based on their hydrodynamic simulations of large- 
scale structures. However, some O VI absorbers show strong 



evidence of hot gas (e.g. Tripp et al. _ 200 ll : I Savage et al. I 



I2OO5I : iNaravanan. Wakker fc Savage ll2009f K Larger statist! 
cal studies employing more robust diagnostics are required 
to reliably interpret the natures of O VI absorbers. 

3C 273 has always been one of the first extragalac- 
tic targets observed with new technology fo r a variety of 
astronomical topics (e.g. lUlrich et al. I Il980l ). In addition 
to 3C 273 being extremely bright, it also lies behind the 
Virgo Cluster, which makes it particularly useful for stud- 
ies of absorber-galaxy relations and the physical condi- 
tions in halo and filament gas, from the cluster out to 
the vicinity of the QSO. In a partial listing of succes- 
sor projects after the early HST papers cited above, the 
sight line was the sub ject of further GHRS observations 
(jWev mann et al. II 19951 ) and a numbe r of galajcy studies re- 
lating to absorber environments (e.g. 'Salpeter fc Hoffman! 
1995; Hoffman ct al. 1998; Grogin, GcUcr, fc Huchra 199^; 



ITmpev. Petrv. fc Flint I 119991 '). ISembach et al. 



(120011) 



plored the Ly^ forest of 3C 273 with FUSE, and lTripp et al. I 
(I2OO2I ) examined the physical states of two high column den- 
sity Virgo Cluster absorbers. iRosenberg et al.~l (|2003l ) com- 
pared Virgo metal systems toward RXJ1230. 8-1-0115, which 
is the nearest bright low redshift QSO (z — 0.117, 55 arcmin 
away), and concluded that a large scale structure filament 
may produce correlated absorption. 

Preliminary results for ~ 7 km s~^ resolution 
STIS echell e spec trosco py of 3C 2 7 3 wer e presented by 
iHeap et al. I (|2002l ') and iHeap et al. I (|2003h . We present a 
full analysis here, including a coverage of Galactic and IGM 
absorption system parameters. We also make use of com- 
plementary data from FUSE, with recent improvements in 
data processing, reduction and analysis. The QSO 3C 273 
lies at Galactic £ — 289.95°, b — 64.36°, and thus serves as 
a very useful probe of the Galactic halo. Over wavelengths 
covered by the Lya forest, our data have S/N — 20 — 30 
per ~ 3.2 km s ~^ pixel with 1.7 pix els per resolution ele- 
ment at 1200 A l|Dressel et al. 1120071 ) . except for small dips 
at the ends of orders. The superior signal and resolution 
of our data permit us to probe the low-z Lya forest to 
log Nhi ~ 12.3 — 12.5, which is the best observation done 
to date, regardless of resolution. This column density limit 
corresponds to ma s s overdensities of log jpH /ph) ~ 0.4 ±0.2 
iDave et al. Ill999l : iDave fc Tripp1l200ll ) as calculated from 
hydrodynamical models in the cold dark matter scenario. 
Similar mass overdensities at 2 ~ 3 correspond to log Nhi ~ 
14.0 — 14.5 (though the true density-column density relation 
is likely complex due to the non-linear growth of large per- 
turbations). The high S/N of our spectrum also allows us to 
detect BLAs with good efficiency. The study of hot halo gas 
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for the Vir go absorbers afong this sight line (|Tripp et al. I 
|2002| .[ 20081 ') prompted us to verify the component structure 
of the z = 0.0053 absorber with the FUSE data. While a 
number of recent Lya forest publications tend to include 
more and more sight lines per paper, our focus on the single 
sight line of 3C 273 merits an in-depth analysis of the Lya 
forest from a special perspective, due to the high quality of 
the spectrum and the value added by existing and future 
deep galaxy surveys in the region. 

This paper is organized as follows. We describe the STIS 
and FUSE observations in § (2] The absorption line selec- 
tion procedure details are in § O as are a list of Galactic 
absorption lines in the STIS data and an upper limit on 
C IV absorption for a high H I column density absorber at 
z = 0.0665. We explain our simulations to constrain the 
probability of feature detection in § 13.41 Our results and 
discussion for the Lya forest, BLAs and their relation to 
O VI absorbers and galaxies, clustering and voids, and the 
z = 0.0053 Virgo absorber velocity structure are in Sjl] Our 
conclusions are summarised in § [5] We adopt a cosmology of 
Ho = 70 km s"i Mpc"\ fib = 0.3 and A = 0.7 throughout 
this work. 



the zero point by cross referencing with HST data for in- 
terstellar lines. Both the FUSE LiF and SiC channel data 
show some improvement in terms of resolution. In the 2001 
reduction, the screening limits for pulse height distribution 
were (manually) set with care, but in this work, the wider 
default values were used. The use of these wider values was 
able to reduce the noise compared to the earlier reduction. 
Aside from a small stretch in the LiF2 data and cases where 
a feature fell on a fibre bundle boundary, the velocity scale 
is generally good to ±5 km s~^, save for a few cases where 
there is an offset of as much as half a resolution element 
(22-25 km s~^). The SiC backgrounds are good except for 
SiCl at A < 920 A. We refer the reader to Sembach et al. I 
for further details. 

For profile-fitting, we originally used night-only FUSE 
data; later we tried using all of the data but found no sig- 
nificant improvement in the results. Continua for the FUSE 
spectra were created using interactive line_norm routines. 
We permitted velocity offsets as a free parameter on the 
order of the FUSE uncertainties as needed to improve indi- 
vidual fits. Absorption system parameters presented in this 
paper use the night-only data. 



2 OBSERVATIONS AND REDUCTIONS 

2.1 STIS spectra 

3C 273 was observed with HST and STIS using guaranteed 
time from Program 8017 to the STIS Instrument Definition 
Team (IDT). The instrumental setup used grating E140M 
for seven orbits (18671 sec) on 2000 May 2 and 2000 June 
21-22, with the 0.2" x 0.2" slit. For the E140M grating, the 
STIS Instrument Handbook gives a dispersion of A/91700 A 
pix~^, a FWHM for the fine spread function for the 0.2" x 
0.2" slit of 1.4 and 1.3 pixels at 1200 and 1500 A respectively, 
and a resolving power for a two-pixel resolution element of 
45800. 

The data were processed and spectra extracted with 
CALSTIS and a suite of programs from the STIS ID Team at 
Goddard Space Flight Center, including corrections for scat- 
tered light and hot pixels and customised merging of echelle 
orders to minimize echelle ripple (written by Don Lindler). 
We constructed a continuum using a combination of routines 
from automated autovp (|Dave. Dubinski. fc Hernguist I 
Il997h and interactive line_norm routines (D. Lindler), be- 
cause the regions around emission lines were best done with 
manual fitting. Longward of our Lya forest sample region, 
two regions at 1402-1405 and 1425-1430 A were lost due 
to vignetting from the STIS repeller wire. The total STIS 
wavelength coverage is 1142-1709 A. Redward of 1616 A, 
there are five small (0.15 — 0.58 A) inter-order gaps. 

2.2 FUSE spectra 

We use FUSE spectra from the principal investigator pro- 
gram P10135. A n analysis of the FUSE data was done by 
ISembach et al. I 12001 ), which employed a customised ver- 
sion of the FUSE pipeline Calfuse 1.8.7. In this study, we 
used similar techniques and Calfuse 3.1.3. We screened the 
data for valid photon events, registered the individual spec- 
tra by cross-correlating on the interstellar features, and set 



3 ABSORPTION LINES 
3.1 Selection 

Absorption line selection was a multi-step process. We be- 
gan the process by pre-selecting an inclusive line list using 
a series of selection criteria based on significance. Our ini- 
tial step was to make a preliminary selection of absorption 
features from the summed STIS data down to the 4.0(7 sig- 
nificance level with a Gaussian filter based on the autovp 
routine, with half- widths of 8, 12, 16, and 20 pixels. The 
large width provides sensitivity to broad Lya absorbers and 
partially resolved complexes. As a further preliminary (and 
complementary) constraint, we then only included signif- 
icant features within the regions fiagged as significant in 
absorption with a simple equivalent width significance cri- 
terion, based on a ~ 3.3(t threshold for contiguous pixels 
below the continuum; such a criterion is more sensitive to 
narrow absorbers. Whenever possible, we confirmed spectral 
features by using the FUSE data. In particular, all (9 out of 
9) absorbers with logNui > 13.24 over 0.003 < z < 0.147 
were confirmed by inspection of Ly/? data, and Ly/3 inspec- 
tions showed no useful information for lower column density 
Lya forest systems. The selection algorithm to this point 
was designed to be largely objective and repeatable, aside 
from continuum fitting around emission features. We then 
visually checked each absorption feature to avoid including 
noise spikes or potential continuum errors. 

We subsequently tested the effectiveness of this selec- 
tion technique by subjecting it to simulated spectra with 
similar resolution to the data to determine the 80% detection 
probability based on a grid of Doppler parameters, H I col- 
umn densities and S/N values. Details are in i|3.3l below. The 
decrease in detection probability drops steeply from nearly 
100% to nearly zero over 0.1 dex in column density for a 
given S/N and Doppler parameter, which enables straight- 
forward determinations of the 80% detection efficiency to 
< 0.05 dex accuracy. Such a steep decline is not unreason- 
able given that the absorbers are of such low column density 
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(log Nhi ~ 12.5) that they are weU on the hnear part of the 
curve of growth. Their equivalent widths change little over 
the limited range of Doppler parameters and S/N (factors 
of 2) which we explored, and their velocity widths also only 
vary over a small interval. Therefore, when detecting weak 
but significant features relying in part on a boxcar approx- 
imation to an absorption profile, a range of 0.3 dex (factor 
of 2) in column density would be a large range to expect a 
uniform detection threshold probability. 

We made a concentrated effort to understand the rate of 
false detections, due to our desire to probe as deeply as possi- 
ble down the H I column density distribution. Therefore, as 
an additional test, we examined the spectrum > 5000 km s~^ 
redward of the QSO's Lya emission all the way to the red 
end of the STIS data (1710 A) for unidentified absorption 
features to test for the frequency of appearance of potential 
spurious features. The most significant one at 1667.9 A had 
a significance 2.3(7, which was well below our preliminary 
selection threshold. It would have Lya fitting parameters of 
log Nhi ~ 12.53±0.13 and Doppler parameter b = 13.4±6.3. 
The relatively large (47%) error in b would make it unlikely 
to be considered a reliable Lya absorber, even if it had some- 
how been included in our preliminary absorber list. The fea- 
ture does not correspond to any plausible metal line for any 
Lya absorber in our sample. The final STIS E140M spec- 
trum of 3C 273, with absorption systems marked, is shown 
in Figure [T] 

3.2 Profile fitting technique 

The data were profile-fitted with vpfitQ dWebb I I 198711 usmg 
Voigt profiles convolved with the local STIS line spread func- 
tion (LSF) taken from the STIS Instrument Handbook. Os- 
cillator strengths are from the vpfit atomic data list (2004 
Dec. version), based on , Morton . (200ii). Multiple transitions 
were used for simultaneous fits whenever they could improve 
constraints on column densities and Doppler parameters. We 
also used, as necessary, the vpfit capabilities to allow for 
offsets in the local continuum level and in velocity for each 
spectral interval used in a fit. The latter would in particu- 
lar compensate for small velocity uncertainties in the FUSE 
data. 



htt p : / / www . ast .cam .ac.uk/^rfc/ vpfit. html 
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Figure 1. 3C 273 data and la errors (red) binned by 3 pixels for 
presentation only. Grey dash- triple dotted curve: 80% detection 
threshold in logAf(H I) (right axis). Long, bold (red) ticks: Ly« 
lines with metals. Long unlabelled (blue) ticks: Lya forest {z on 
upper axis) with log Nhi < 12.5 in the complete survey (? [2)l. 
Thick (blue) tick: reliable BLA. Short (blue) ticks: other Ly« 
lines and higher order Lyman lines (labelled). Heavy dotted (red) 
ticks: intervening metal lines (labelled). Medium dotted (green) 
ticks: Galactic metal absorption. 
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Figure 1. Continued. 
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3.3 Simulated spectra to determine detection 
probabilities 

To give us a clearer picture of our Lya line detection prob- 
ability as a function of S/N ratio and Doppler parameter, 
we analysed 880 simulated Lya lines. The emphasis here 
is on characterizing the detection and recovery of of weak 
absorbers and their line parameters, because the 3C 273 
data uniquely probe the low end of the Lya forest col- 
umn density distribution. The simulations were in addition 
to the ones done for PKS 0405-123 described in Paper I. 
The simulated lines were generated using the STIS LSF 
and a grid of values for S/N, Doppler parameter and H I 
column density. The S/N values were set to either 20 or 
40 per pixel, Doppler parameters were set to ranges with 
averages of 17.5, 22.2 or 24.7 and 35.0 km s~^, and the 
input log Nhi values varied over 12.21-12.80 for S/N=20 
and 11.85-12.55 for S/N = 40. The simulated spectra were 
continuum-normalised, so any errors, systematic or random, 
involving continuum fitting are not reflected in the subse- 
quent simulation analysis. There were also no noise spikes 
inserted into the spectra. We calculated statistics comparing 
the line parameters for simulated lines which were success- 
fully recovered vs. the input values. Results for the vari- 
ous combinations of input S/N ratio and Doppler parame- 
ter are shown in Tables [T] and [H The mean recovered mea- 
sured Doppler parameters and column densities match the 
mean input values within the mean of the profile fitting 
errors. The marginal tendency for the narrowest lines to 
be recovered with higher Doppler parameters is likely due 
to the effects of noise broadening input for features near 
the detection limit for recovered lines; noise which would 
make lines narrower presumably would cause the features 
to drop below the (essentially equival ent width-defined) de- 
tection threshold (|Rauch et al. lTl993l ). The recovery of the 
line parameters (e.g., the 6- value) was found to be even bet- 
ter for stronger lines detected at high significance, up until 
the point where the H I lines start to saturate badly. The 
spurious detection rate is < 1% among all the simulations 
above the 80% detection probability threshold, and thus is 
not considered significant. The boundary of the 80% detec- 
tion probability threshold can be parametrically fitted by 
logNni = 12.11 + 0.271og(6/24.7)-1.121og(snr/40.), where 
snr is the signal to noise ratio (S/N) per pixel, and estab- 
lished using simulated lines with 16<6<36 km s~^. The 
process was adapted from the method used in Paper I. 



We consider this parametrization useful to determine our 
80% detection threshold over a Doppler parameter range of 
10-40 km s~^, with the bounds based both on our simu- 
lations and to compare with the analysis of iLehner et al. I 
(|2007al ). 



2 The coefficient of 0.27 for Iog(6/24.7) for this paper is ~ 25% 
lower than the analogous value of 0.34 from Paper I, and the 
snr coefficient in this work is 11% higher. The difference in the 
log(6/24.7) coefficient may be caused by the relatively lower col- 
umn density range used for the input lines for this study, in that 
the limited Doppler parameter range used affects the equivalent 
widths and thus detectibility of the lines less strongly. The differ- 
ence in the snr coefficient may also stem from a similar cause, or 
from either the higher S/N in the 3C 273 data and/or the use of 
only two input S/N values in this study, 20 and 40. In comparison, 
for PKS 0405-123, we used three evenly distributed in log space 
over the range we tested. The difference in results between a coef- 
ficient of 1.12 and 1.00 (without adjusting the other coefficients) 
is only 0.03 dex for our minimum S/N of 21, and the current 
value of 1.12 leads to a more conservative detection threshold as 



snr/40 < 1. We do not believe that the results of this paper are 
significantly affected by the difference in parametrization. 

We will publish an erratum for Paper I (Williger et al. , in prep.) . 
In that work, a subset of absorbers principally below the 80% 
detection threshold had erroneously high calculated significance 
because the contribution of continuum uncertainties to the total 
equivalent width errors was not included in the software. 
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Table 1. Statistics of simulated data: Doppler parameters 



S/N" Doppler parameter (km s ^) Doppler parameter (km s ^) Doppler parameter (km s ^) 

input*" recovered'^ input'' recovered'^ input'' recovered'^ 

20 17.2 <6< 17.8 19.0 ±5.0(4.1) 22.2 ±5.5 21.8 ± 8.0(5.5) 35.0 ± 0.6 35.3 ± 8.1(8.1) 

40 17.2 <6< 17.8 19.5 ± 6.8(4.5) 24.7 ±0.1 26.7 ± 7.0(6.3) 35.0 ± 0.6 35.0 ± 7.0(7.2) 



" Signal to noise ratio per pixel in simulated STIS spectra. 

Mean and first order about the mean for input Doppler param- 
eters in simulations which were recovered. The simulations for 
S/N=20, (6) = 22.2 km s~^, were made from a distribution made 

at the beginning of our analysis, with a much wider standard 
deviation (5.5 km s^-'^) in input b than the others, hence the 
difference with the (6) = 24.4 km s^-"- case (the geometric mean 
of 17.5 and 35.0). The results between the two cases are within 
the scatter, so we did not maJce a simulation run with S/N=20, 
(b) = 24.7 km s"!. 

° Mean and first order about the mean for recovered Doppler 
parameters in simulations, with the mean in the formal Icr 
profile fitting errors in parentheses. The mean recovered Doppler 
parameters are all within the mean profile fitting errors of the 
input values. 



Table 2. Statistics of simulated data: H I column densities 





Doppler'' 


log 


;iV(H I) 




parameter (km s"'^) 


input° 


recovered'' 


20 


17.5 


12.54 ±0.12 


12.55 ±0.12(0.07) 


40 


17.5 


12.17 ±0.12 


12.20 ±0.11(0.09) 


20 


22.2 


12.41 ±0.10 


12.41 ±0.13(0.08) 


40 


24.7 


12.26 ±0.12 


12.28 ±0.13(0.08) 


20 


35.0 


12.64 ±0.12 


12.65 ±0.13(0.08) 


40 


35.0 


12.36 ±0.12 


12.37 ±0.13(0.07) 



" Signal to noise ratio per pixel in simulated STIS spectrum. 

Simulation sets are specified by input Doppler parameter. 
° Mean and first order about the mean for input log H I column 
densities in simulations which were recovered. 

Mean and first order about the mean for recovered log H I 
column densities in simulations, with the mean in the formal Icr 
profile fitting errors in parentheses. 
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3.4 Profile fits to data 

In addition to H I lines, we find and list, for completeness, 
profile fits or column density limi ts (from vpfit or from the 
apparent optical depth method. ISavage &: Sembach Ill99ll ) 
for a number of Galactic and intervening metal absorber 
species. 

Galactic absorption. We find Galactic absorption from 
H I, C I, G r, C II, G ir, G IV; N I, N V, O I, Mg II, 
Al II, Si II, Si III, Si IV, P II, S II, S III, Mnll, Fe II, Nill. 
Galactic absorption ranges over —135 < t; < 68 km s~^, con- 
sistent with t he high velocity absorp tion found in the FUSE 
spectrum bv lSembach et alTI (|200ll '). A complete list of ab- 
sorber profile fitting parameters is in Table [S] In the case 
of saturated lines, we calculate lower limits by th e apparent 
optical depth method . Savage fc Sembach l (Il99ll'). whi ch are 



consistent with the results of ISembach et al. ( 200 ll ). The 



maximum optical depth used is a function of local S/N and 
continuum error {rmax ~ ln[\/2/(S/A'^)] per pixel), in which 
both effects contribute to similar degrees. For the 3G 273 
data, Tmax ^ 2.3 3.0. 

For two of the highest Nhi column density systems 
toward 3G 273, members of this collaboration and others 
have investigated metallicit y limits and und e rtaken stud- 



ies in the low z IGM (e . g. loave et al. |[200ll : iTripp et al~l 



l2002l : IStocke et al~ll2007l : iTripp et al. II2008I ). Such low red 



shift, high column density absorbers are extremely useful for 
comparison with nearby low luminosity faint galaxies e.g. 
the dwarf galaxy found 71 kpc away from the z = 0.0053 
Virgo cluster absorber (|Stocke et ain[2"oo4 ). For complete- 
ness and comparison with their results, we provide here up- 
per limits for C IV for each of these Nhi absorbers, and for 
a sum of the C IV regions over a consistent —100 < Av < 
100 km s~^, using the above-mentioned apparent optical 
depth metho d. The z = 0.003369 absorber was found by 
iTripp et al. I to have rest equivalent width Wr,civ < 27 mA 
over -100 < Av < 100 km s-\ or logiV(C IV) < 12.8 
assumin g an absorber on the linear part of the curve of 
growth. IStocke et all found log iV(G IV) < 12.27 (3cr) for 
the z = 0.066548 system over -50 < Av < 50 km s"\ 
We measured an upper limit for the z — 0.066548 absorber 
of Wr,civ<50 mA (3a) over -100 < Av < 100 km s-\ 
or logiV(G IV) < 13.10 We stacked the data for the two 
systems and found MKr,C7V < 19 niA (3(t) -100 < Aw < 
100 km s"\ or logAr(G IV) < 12.70 



detection threshold for our worst case S/N = 20 in the Lya 
forest and assuming a Doppler parameter of 6 = 40 km s~^ 
0. This low H I column density threshold sample covers the 
range 0.020 < z < 0.139 and contains 21 Lya absorbers. We 
excluded absorbers within Av ~ 5000 km s~^ of 3G 273 to 
minimize the impact of the proximity effe ct ([Murdoch et al. I 
ll986l : lBaitlik. Duncan fc Ostriker Ill988h . 

In later subsections, we will discuss in detail our exam- 
inations of the Lya forest Doppler parameter distribution, 
redshift density dN/dz, broad line population and column 
density distribution. We do not distinguish between metal 
and Lya-only systems in this analysis, the same approach 
as in Paper I. Metals have been found in progressively lower 



column density Lya systems at z > 2 (e. g. Gowie et al. Il99^ : 
ISongaila fc Gowie lll99HPettini "etlj~ll200ll . and references 
therein). Also, the absorbers we can detect at low redshift 
probably a re best related to large density perturbations at 
z>2 (e.g. ISchave I [ 200 ll ) which are commonly associated 
with metals. Finally, the limited wavelength range of the 
data does not permit a uniform coverage for metal detec- 
tion (in particular for G IV, limited to 2:<0.1). 



4 RESULTS AND DISCUSSION 
4.1 Lya absorber statistics 

We use our detection probability results to define a sample 
for log(iV(II I)) > 12.5, which is approximately the 80% 

* The difference in upper limits arises in that the IStocke et al. I 
limit is for one resolution element, whereas our limit is over 
200 km s~^ in veloci ty space or ~ 28 resolution elements, to 
be consistent with the lTripp et al.1 measurement. 
^ This limit is slightly tighter than the 24 mA result expected 
from stacking the spectra according to the inverse square of the 
S/N. The difference in our result may be from the way each team 
determined the local continuum and S/N, in particular as these 
were measured over two different velocity intervals. 



° Note that by choosing b = 40 km s~ , this leads to a higher, 
more conservative 80% probability detection threshold in Nhi 
by 0.06 dex than by choosing 6 = 25 km s~^, which is typical 
of high resolution, high S/N Lyo forest data. We prefer to be 
conservative, given the possibility of unresolved blends caused by 
our limits in the achievable S/N. This choice of b does not affect 
our subsequent analysis in any way except for this small change 
in the threshold for log Nhi. 
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Table 3. Galactic absorption lines toward 3C 273 from STIS data 



species 


z 


error" 


(km s ) 


error" 


log A'^ (km s ^ ) 


error" 


Imes (A) 


O I 


-0.000451 


0.000041 


61.5 


22.1 


13.55 


0.12 


1302 


C IV 


-0.000240 


0.000003 


11.6 


1.8 


13.45 


0.09 


1548,1550 


Si IV 


-0.000230 


0.000003 


7.9 


1.4 


12.57 


0.06 


1393,1402 


Si II 


-0.000224 








> 13.24 




1190,1193,1260,1304,1526 


O I 


-0.000220 








> 13.59 




1302 


S II 


-0.000159 


0.000016 


10.1 


4.9 


13.74 


0.24 


1250,1253,1259 


Si IV 


-0.000103 








> 13.19 




1393,1402 


S III 


-0.000084 








> 14.26 




1190 


Fe II 


-0.000071 








> 14.71 




1608 


N I 


-0.000067 








> 15.04 




1199,AA1200 


Mn II 


-0.000067 




14.2 




12.60 


0.19 


1197,1199, 1201^* 


Ni II 


-0.000067 


0.000004 


22.3 


1.7 


13.73 


0.03 


1317,1370,1454 


C II 


-0.000067 








> 15.54 




1334 


C IV 


-0.000066 


0.000003 






> 14.46 




1548,1550 


S II 


-0.000060 








> 15.22 




1250,1253,1259 


Si IV 


-0.000060 


0.000012 


48.8 


2.3 


13.73 


0.04 


1393,1402 


P II 


-0.000056 


0.000017 


29.1 


7.5 


13.69 


0.10 


1152,1532 


Si II 


-0.000067 








> 15.03 




1190,1193,1260,1304,1526 


C II* 


-0.000045 


0.000002 


14.3 


0.8 


13.83 


0.02 


1335 


C I* 


-0.000044 


0.000000 


9.1 


0.0 


12.24 


0.26 


1561,1656,1657'= 


C I 


-0.000044 


0.000005 


9.1 


2.4 


12.79 


0.07 


1277,1280,1328,1560,1656 


Si III 


-0.000042 








> 14.02 




1206 


Mg II 


-0.000038 


0.000012 


27.9 


5.1 


15.43 


0.06 


1239,1240 


N V 


-0.000037 


0.000005 


59.8 


2.2 


13.93 


0.01 


1238,1242 


Al II 


-0.000033 








> 13.62 




1670 


H I 


-0.000027 


0.000005 






20.23 


0.00 


Lyo,/9 


O I 


-0.000025 








> 15.38 




1302 


Fe II 


0.000068 








> 14.28 




1608 


N I 


0.000068 








> 14.66 




1199,AA1200 


Mn II 


0.000068 




4.4 




12.63 


0.18 


1199,AA1200'* 


S II 


0.000069 


0.000001 


9.4 


0.6 


14.75 


0.01 


1250,1253,1259 


Ni II 


0.000074 


0.000004 


13.3 


2.0 


13.41 


0.05 


1317,1370 


C II* 


0.000078 


0.000001 


7.6 


0.7 


13.73 


0.02 


1335 


C I 


0.000081 


0.000001 


6.4 


0.7 


13.27 


0.02 


1277,1280,1328,1560,1656 


C I* 


0.000081 


0.000000 


6.4 


0.0 


12.67 


0.09 


1561,1656,1657= 


C IV 


0.000138 


0.000012 


20.7 


4.4 


13.32 


0.15 


1548,1550 


Si II 


0.000159 








> 14.59 




1190,1193,1260,1304,1526 


S II 


0.000171 


0.000009 


1.3 


4.1 


13.13 


0.21 


1250,1253,1259 


O I 


0.000190 


0.000015 


13.0 


5.7 


13.41 


0.19 


1302 


C II 


0.000206 








> 14.89 




1334 


Si IV 


0.000225 


0.000020 


15.2 


11.6 


12.21 


0.56 


1393 



" Errors are 1 a, and assume that the component structure is 
correct. 

^ Upper limits from VPFIT or the apparent optical depth method 
of ISavaee fc Sembach I lll99l[) . 
Lines used in profile fits. 

Mn II tied in redshift and Doppler parameter to N I. 
''CI* tied in redshift and Doppler parameter to C I. 
f We find an unidentified feature at 1149.3 A (also in the FUSE 
data) . 

f We find a weak feature at 1304.74 A to be an artefact just 
redward of an echelle overlap region, and not securely identifiable 
with any plausible ISM or IGM absorption. 
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4-1.1 Redshift density and intervening absorber list 

Measuring the Lya forest redshift density has been one of 
the biggest goals among HST Key Projects, and initially 
yielded surprising re sults. From among t he first resul t s with 
FOS Bahcall ct ah] l|l99ll ) and GHRS Morris et al. I ([l99l|) 
toward 3C 273, the Lya redshift density was found to be 
higher than expected from extrapolations of ground-based 
data. To compare results with recent high resolution data 
in the literature, we present data for 13.1 < log Nhi < 
14.0 and 14.0 < log Nhi < 17.0, using the restrictions of 
6 < 40 km s"'^, and errors in b and Nhi of < 40% for 
both 3C 273 and the l i teratu re values (to use the same cri- 
teria as iLehner et al. I l2007al ) . We also present results for 
log Nhi > 12.5. Due to the low redshift of 3C 273, the data 
are unaffected by losses from inter-order gaps in the STIS 
data. Similar to our analysis in Paper I, we excluded pix- 
els from metal absorption with optical depth r > 2, which 
removed Az = 0.002 from the samples for log Nhi > 13.1 
and log Nhi > U.O^ This left an unblocked redshift path 
of Az = 0.135 for log Nhi > 13.1 and for log Nhi > 14.0. 
In an analogous way, we find Az — 0.117 for log Nhi > 12.5 
(which requires a smaller, higher S/N subset of the data). 
The absorption di stance AX = 0.124 (iTvtler IITqStI ). using 
the prescription of lLehner et al. I l|2007a ). which will be used 
for comparing H I column density distributions in the next 
section. 

We list in Table |4] intervening If I and related metal 
redshifts, column densities and Doppler parameters, includ- 
ing several systems either outside of the unblocked redshift 
intervals corresponding to the sensitivities mentioned above, 
or at 12.3 < log Nhi < 12.5 which nevertheless fulfill all of 
the other selection criteria. We show in Figure[2]the Doppler 
parameters and H I and column densities for all Lya forest 
systems that we detected using the procedure described in 
§ |3l overlaid with the our parametrically fitted 80% detec- 
tion sensitivities. 

4.1.1.1 A. High column density lines log Nhi > 14.0 

For absorbers with log Nhi > 14.0 from among our full 
sample of 21 systems, we find four over the interval 0.002 < 
z < 0.139 {logdN/dz = 1.471o;3q). ^'^^ restrict our sample 
to 6 < 40 km s~^ and errors in b and Nhi to < 40%, the 
number reduces to 3 (log dN/dz = 1.351q gg). The difference 
arises from the lower Nhi component of the z = 0.0053 
Virgo absorber, which has a column density error of a 
factor of 2 and is examined in more detail in § 14.41 We 
compare our results with the literature (Fig. |3} , and illus- 
trate the scatter for seven other sight line s with contigu- 
ous redshift coverage of Az > 0.1 3 from IWiUiger et alT] 
l|2006l . PKS 0405-123, Az = 0.413)Fl [Sembach et al. I l|2004l . 
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0.020<z<0.139 
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log N(HI) 

Figure 2. H I column densities and Doppler parameters. Lya 
forest sample over 0.020 < z < 0.139, which is the redshift range 
used for our sample of log Nhi ^ 12.5. Open diamonds: data 
points. For clarity, only error bars for every second point are 
shown. The line indicates the parametric approximation to the 
80% detection probability from simulations in § 13.41 based on 
S/N> 21. Three absorbers with 12.3 < log Nhi < 12.5 used for 
column density studies are also included. 



PKS11164-215, Az = 0.1 33) and iR ichter et aTl (|2004l . 
PG 12594-593, Az = 0.247). iLehner et a l. (20071hE0226- 
4110, H1821-h643 ar id PG0953-F4 1 5, Az = 0.397, 0.238, 
0.202 respectively), lAracil et al.1 (|2006l . HS 0624-^6907, 
Az = 0.329). The redshift density is marginally high (sim- 
ilar to PKS 0405-123) but still consistent with other mea- 
surements from the literature shown in Fig. [3] (right panel), 
similarly selected for & < 40 km s~'^and with errors in b and 
Nhi of < 40%. The small number statistics and large cos- 
mic scatter (~ 0.7 dex in dN/dz for individual sight lines) 
make make further analysis difficult without larger samples. 
The cosmic scatter range appears to persist to z ~ 1. To 
illustrate the effect of the restriction in b and in the er- 
rors, we also included the redshift density f r om th e large, 
lower resolution sample of IWevmami et al. I j 19981 '). which 
was selected on the basis of rest equivalent width and can- 
not filter out & > 40 km s~^ systems, and is ~ la higher 
than our result. We note th a t the z ~ dN/dz determi- 
nation of IWakker fc Savage ] (|2009D is consistent with the 
IWevmann et al. result, as both studies relied on equivalent 
widths rather than column densities with restrictions on the 
Doppler parameters. The redshift densities for restricted, 
high resolution 2<0.3 sight lines are in general a factor of 
--^2 — 3 lower than the analogous 1 < z < 2 sight lines. 



^ Metal absorption from pixels with r < 2 can also reduce the 
probability of detection an IGM Lya absorber, but less so than 
those with t > 2. The value of r = 2 is a compromise to elimi- 
nate the worst metal absorption regions without losing too many 
pixels with recoverable Lya data. Other factors such as multiple 
transitions from a metal ion can mitigate the masking effect of 
metal absorption. For example, we recovered the presence of Lya 
absorption blended with Galactic Si IV. 

* We corrected the line list for spurious absorbers, and will pub- 
lish the updated list in an erratum (Williger et al. 2009, in prep.). 
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Table 4. Intervening absorption lines toward 3C 273 from STIS 
data 



species 


z 


error'^ 


(km s ) 


error" 


log N 


error" 


imes (A) ' 


H I 


0.002630 


0.000039 


125.8 


44.9 


13.04 


0, 


.14 




H I 


0.003369 


0.000002 


37.4 


0.8 


14.22 


0, 


.02 


Lyo, 


H I 


0.005251 


0.000021 


26.1 


1.3 


14.33 


0, 


.31 


Lya, ^= 


Si II 


0.005266 


0.000005 


5.2 


2.6 


11.78 


0, 


.09 


1190,1260 


H I 


0.005295 




15.0 


0.0 


15.68 


0, 


.04 


Lya, B, ■7, (5, e, cT, 0, 


C II 


0.005295 




8.6 


0.1 


12.80 


0, 


.07 


1036,1334 


Si III 


0.005295 


0.000002 


8.2 


2.1 


12.36 


0, 


.07 


1206, sets z for Si III, H I 


H I 


0.007165 


0.000029 


38.8 


12.1 


12.67 


0, 


.12 


Lya 


H I 


0.007588 


0.000037 


64.0 


18.0 


12.86 


0, 


.10 


Lya 


H I 


0.026225 


0.000006 


26.7 


2.5 


12.95 


0, 


.03 


Lya 


H I 


0.029348 


0.000028 


17.1 


6.4 


12.77 


0, 


.28 


Lya 


H I 


0.029471 


0.000008 


19.7 


7.8 


13.24 


0, 


.16 


Lya'^''= 


H I 


0.029578 


0.000016 


11.6 


5.9 


12.63 


0, 


.30 


Lya 


H I 


0.032798 


0.000006 


22.8 


2.6 


12.82 


0, 


.04 


Lya 


H I 


0.039373 


0.000012 


28.2 


5.2 


12.63 


0, 


.06 


Lya 


H I 


0.046639 


0.000013 


24.4 


5.5 


12.51 


0, 


.08 


Lya 


H I 


0.049000 


0.000002 


28.7 


0.7 


13.56 


0, 


.01 


Lya'^''= 


H I 


0.049953 


0.000004 


16.2 


1.7 


12.83 


0, 


.03 


Lya 


H I 


0.054399 


0.000017 


23.6 


6.6 


12.41 


0, 


.10 


Lya 


H I 


0.063510 


0.000018 


17.9 


3.9 


12.96 


0, 


.14 


Lya 


H I 


0.064078 


0.000025 


39.8 


10.9 


12.62 


0, 


.09 


Lya 


H I 


0.066548 


0.000001 


36.7 


0.5 


14.08 


0, 


.01 


Lya, /3<= 


H I 


0.073738 


0.000034 


64.0 


15.0 


12.74 


0, 


.08 


Lya 


H I 


0.074286 


0.000033 


51.2 


13.8 


12.62 


0, 


.10 


Lya 


H I 


0.075377 


0.000009 


27.1 


3.8 


12.77 


U 


.Uo 


Lya 


H I 


0.085882 


0.000028 


34.9 


11.3 


12.41 


0, 


.11 


Lya 


H I 


0.087632 


0.000008 


43.3 


3.1 


13.09 


0, 


.03 


Lya 


H I 


0.089746 


0.000012 


30.5 


4.4 


13.08 


0, 


.06 


Lya 


H I 


0.089898 


0.000009 


12.8 


4.3 


12.60 


0, 


.15 


Lya 


H I 


0.090110 


0.000004 


30.1 


1.8 


13.29 


0, 


.02 


Lya'*'= 


H I 


0.109141 


0.000015 


18.8 


6.2 


12.33 


0, 


.11 


Lya 


H I 


0.109380 


0.000009 


18.7 


3.8 


12.56 


0, 


.06 


Lya 


O VI 


0.12003 








13.45 


0, 


.10 


1032,1038-'' 


H I 


0.120041 


0.000002 


21.5 


0.7 


13.50 


0, 


.01 


Lya'^'S 


H I 


0.141676 


0.000024 


32.6 


9.6 


12.36 


0, 


.10 


Lya 


H I 


0.146575 


0.000005 


39.0 


1.7 


14.08 


0, 


.03 


Lya, 


H I 


0.150204 


0.000021 


32.9 


7.8 


12.37 


0, 


.08 


Lya 


H I 


0.157787 


0.000006 


21.5 


2.1 


12.63 


0, 


.03 


Lya 


O VI 


0.15779 








13.39 


0, 


.11 


1032-f 



" Errors are 1 cr, and assume that the component structure is 
correct. 

^ Lines used in profile fits. 

We find an unidentified feature at 1149.3 A (also in the FUSE 
data). 

Inclusion of Ly/3 makes no significant diff'erence to the fit. 
Lya absorption con firmed by Ly/3. 
■^ iTripp et all | |2008^ . 
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Figure 3. Left: Lya forest redshift density dN/dz for 13.1 < 
log Nh I < 14.0. Right: 14.0 < log Nhi < 17.0. Error bars show 
redshift bins and Poissonian errors in dN/dz. We restrict the data 
to b < 40 km s^^, and to errors in b and Nhi of < 40% for 
3C 273 and the l iterat ure values in both cases except for the 
IWevmann et al. I l|l998i ) values, which are based on equivalent 
widths and are shown for comparison. For the open crosses of 
Ijanknecht et al. I l|2006f l. dark symbols denote the highest qual- 
ity data and l ight (re d) symbols all data, as defined in § 4.2 of 
iLehner et aTl l|2007al 'l. 
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4.1.1.2 B. Medium column density lines 13.1 < 
log Nhi < 14.0 In both in our full sample and our re- 
stricted sample for & < 40 km s~^ and errors in b and 
Nhi of < 40%, we find 4 absorbers over the interval 
0.002 <z< 0.139 with 13.1 < log Nhi < 14.0 {log dN/dz = 
1.471q'3q). compare our results with the literature, and 
illustrate the scatter for seven other sight lines with con- 
tiguou s redshift coverage of Az > 0. 13 from Williger et al.] 



l|2006l. PKS 0405-123, Az = 0.212). ISemb"ach et al. I 
PKS1116-h215, Az = 0.133) and iRichter et al 



2004 



2004 



PG 1259-F593, Az = 0.247'). iLehner et aL I (|2007al . HE0226- 



4110, H1821-H643 ar id PG0953-H4 1 5, Az = 0.294, 0.238, 
0.202 respectively), lAracil et alTI (|2006|. HS 06244-6907, 
Az — 0.329). We note that Lehner et al. I estimated com- 
pleteness to log Nhi = 13.2, so the line densities for the 
three objects from that study should be a lower limit. For 
PKS 0405-123, we use a revised line list for PKS 0405-123 
(Williger et al., in prep.) and the redshift interval known 
to be complete to log Nhi = 13.1 (Paper 1). We consulted 
with N. Lehner about the signal-to-noise characteristics of 
the other sight lines in their work (private communication, 
2007). We have split the HS 0624+6907 sight line into high 
and low z halves to keep the redshift intervals to similar 
length compared to that of 3C 273. Results are plotted in 
Fig.|l(left panel). 

The cosmic scatter covers about 0.2 dex in dN/dz at 
2:<0.3, for all sight lines except 3C 273, which is ~ 1.5(t 
below the mean of the other seven sig ht lines. For com- 
parison, the z ~ redshift density of IWakker fc Savage] 
of log dN/dz ~ 2 ± 0.1 (based on their fig. 5b) is just con- 
sistent with the maximum of the individual sight line val- 
ues, again possibly showing the difference in redshift den- 
sity using equivalent widths vs. column densities with re- 
stricted Doppler param eters. At z ~ 1, a larger sample 
jjanknecht et al.HbOOd ) indicates a cosmic scatter range of 
0.35 dex which spans both the 3C 273 value and the other 
seven low z sight lines, so the 3C 273 sight line may simply 
represent the lower end of the redshift density range. The 
high signal and high resolution in the data here make it un- 
likely that blending and noise effects are mainly responsible 
for the paucity of Lya forest lines toward 3C 273 in this 
column density range. The Cosmic Origins Spectrograph 
should reveal more definitively whether logdN/dz ~ 1.5 is 
representative of one of the sparsest parts of the Lya forest, 
or if even more poorly populated sight lines exist. 

4.1.1.3 C. Low column density lines log Nhi > 12.5 

Given the very high signal-to-noise ratio in the 3C 273 data, 
we also calculate the redshift density for log Nhi > 12.5. 
We find 21 absorbers above this threshold over 0.020 < 
z < 0.139, for logdN/dz = 2.25l[J;?|^. To compare the 
3C 273 sample to the m uch larger observational data set 
of (jLehner et al. Il2007al ). we consider separately those sys- 
tems with errors in Doppler parameter and Nhi under 
40% (based on the lower limit la errors for Nhi) and 
6 < 40 km s~^. These criteria focus the sample on rela- 
tive ly cool IGM ga s by excluding broad L ya systems (BLAs, 
e.g. iRichter et al.l i2006a : Richter et al. I 2006b, ') , which are 
interpreted as the signatures of physically distinct low-z 
hot IGM gas. The restricted 3G 273 sample is 18 systems 
(log dN/dz 

density with the full sample. 



2.191qj2), which consistent to la in redshift 



There are few comparisons down to a threshold of 
log Nhi > 12.5 f o r abs orbers at low z. Observationally, 
iDanforth fc ShuU I (|2008| ) found logdN/dz ~ 1.5lo.i for 
z < 0.4 (see their fig. 3), wh ich is 2 — 3g below our v alue 
of logdN/dz = 2.19lo:}2, and lWakker fc Savage I (|2009l ) de- 
termined logdN/dz ~ 2.4 ± 0.2 for z ~ 0, which is con- 
sistent with our value. These results could be understood 
in that the weakest lines in the sample dominate the num- 
ber counts, based on the column density distribution, and 
should affect an e quivalent width based sample (such as 
IWakker fc SavageT) less than a sample with a higher col- 
um n density floor o f log N hi ~ 13.1. From numerical mod- 
els, |Paschoseiar] liroi) calculated logdN/dz ^2.1 for 
the same H I threshold at z — 0.1, which is consistent with 
our measurements to 1 — 1.5a. Less directly, at z — 0, our 
H I column density thres hold corresponds to over densities of 
logipn/pn) ~ 0.4 ± 0.2 (jPave et ainil99^ : IPav c & Trim] 
l2001f ) as calculated from hydrodynamical models in the 
cold dark matter scenario. At z ~ 3, the same overdensity 
level corresponds to log Nhi ~ 14.0 — 14.5, a rough esti- 
mate for which can be read from the right panel of Fig. |31 
dN/dz ~ 2.2± ~ 0.1. Our results are therefore consistent 
with overdensities ollog{pH / pn) ~ 0.4±0.2 within 1.1 cr. We 
caution that the correspondence of overdensity and dN/dz 
may only result from the considered density perturbations 
being weak, as non-linear growth would be expected over 
time to skew the overdensity distribution to high values; the 
low z IGM is relatively complex, containing a significant 
fraction of shock-heated gas (which is related to the BLA 
population) . 

The observation of many more sight lines with COS 
would likely reveal better the extent of cosmic scatter, and 
provide further compari sons between the L ya forest and 
galaxy environments e.g. IChen et al. I (|2005| ). In particular, 
COS could provide a large increase in sight lines probed to 
log Nhi ~ 12.5, though at lower resolution than this study. 

4.1.2 Doppler parameter distribution 

The Doppler parameter distribution for the entire sample of 
21 absorbers has mean, median and standard deviation of 
28, 27, 13 km s~^, respectively (Fig. 3}. In the same man- 
ner as Paper I, we tested whether weak line blending is the 
cause. We divided the strong sample into high and low H I 
column density halves, which occurs at log Nhi = 12.80, 
and performed a Kolmogorov-Smirnov test to determine 
the likelihood that the Doppler parameters are drawn from 
the same distribution. The probability that the samples are 
drawn from the same distribution is P = 0.96. The clos- 
est observed Lya forest absorber pair in velocity space is 
A?; = 31 km s-^ (between z = 0.029471,0.029578), or 4-5 
times the resolution of the STIS data. The next closest pair 
is Aw = 41 km s"\ This can be compared to the mean ve- 
locity interval expected between absorbers (in the absence 
of clustering, which is expected to be weak and shows no ev- 
idence to the contrary, §113]) of {Av) ~ 1500-1700 km s"S 
based on our value of dN/dz for log Nhi > 12.5 over 
0.020 < 2 < 0.139. We would then expect to first approxima- 
tion that there is a ~0.5% chance of an absorber randomly 
being ~ 8 — 9 km s~^ from another. At such a velocity differ- 
ence, the lines could potentially be unresolved, depending on 
the local S/N and the line column densities and Doppler pa- 
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Figure 4. Lya Doppler parameter distribution for a series of 
minimum log Nhi thresholds in the complete sample. 



rameters, with the greatest possibility of blending for broad, 
shallow lines in relatively noisy parts of the data. For the 
high S/N in the 3C 273 data and the & < 40 km s"^ ab- 
sorbers with < 40% errors upon which we are focused, we 
infer that line blending is not significant for the Lya forest 
toward 3C 273. 

We compared re s ults w ith the much larger low z sam- 
ple of iLehner et al. I (|2007al ). They restricted their sample 
to & < 40 km s~^, to correct for the non-Gaussianity of 
the Doppler parameter distribution. If we make the same 
restriction for our data, the sample would have 18 ab- 
sorbers, with m ean, median a n d stan dard deviation 24, 24 
and 8 km s~ \ iLehner et akl l|2007ah found values of 27, 
27 and 8 km s~^, which is c onsistent with our results. The 
log > 13.7, 2 < 0.23 sample of lShull et all (|2000l 'l is also con- 
sistent (mean, median, standard deviation 31, 28, 7 km s~^ 
respectively) . Finally, our b < 40 km s~^ sample is in accord 
with the Doppler parameter mean ofH ave fc Tripp I (|200ll ') 
(25 km s~^), which is complete for log Nhi > 13.0 and in- 
cludes absorbers down to log Nhi ~ 12.6, and thus shows 
agreement wit h their comparison model of a A-dominated 



CDM universe llDave et al. II 19991). ar id aga in with the more 



recent simulations of Paschos et al. I |200^ 



(|Lehner et al. |[2007al ). which contains 270 Lya systems at 
z < 0.44, covering absorption distance Az — 2.40, and a 
stated completeness to log Nhi = 13.2. They calculated the 
differential column density distribution for a subsample with 
fe < 40 km s^^, and errors in b and Nhi to < 40%, for a 
total of 132 absorbers. If we select from our sample in a 
way that is consistent with these rules, the number of in- 
cluded 3C 273 absorbers at 12.5 < log Nhi < 13.2 is re- 
duced to a subsample of 11. We ca lculate a maxim um likeli- 
hood fit of /3 = 1.79 ± 0.07 for the ILehner et alH subsample 
over 13.2 < log Nhi < 16.5, with the distribution consistent 
with a power law to ~ 2a (Kolmogorov- Smirnov proba bility 
P = 0.066). This is consistent with the ILehner et al. I value 
of P — 1.76 ± 0.06. We wish to know wh ether the weak est 
absorbers below the completeness limit of ILehner et al. I are 
consistent with a single power law distribution. The 3C 273 
sam ple would com prise only 7.6% of a combined sample with 
the ILehner et al. I data, over only 4.9% of the absorption 
distance, so it would not contribute much weight to a com- 
bined sample. We theref ore simply extr apolate our best sin- 
gle power law fit for the ILehner et al. I subsample, and plot 
two data points from the complete 3C 273 sample cover- 
ing 12.5 < log Nhi < 13.2 in Fig. O The largest deviation 
shown by the 3C 273 data is in the 12.5 < log Nhi < 12.8 
bin, which is within 0.8cr of the extrapolated fit. The low col- 
umn density end of the column density distribution is thus 
consistent with that at 13.2 < log Nhi < 16.5. 

We can probe to 0.2 dex lower column density by cor- 
recting for the probability of detection at 12.3 < log Nhi < 
12.5 of ^ 50%. U sing the same selection criteria as 
ILehner et all (|2007al ). there are 3 absorbers at 0.020 < z < 
0.139 toward 3C 273 in this extended column density range. 
We plot both the observed point, which shows signs of a 
turnover in dN/dNni , and a factor of two correction, which 
is still consistent to ~ 1.5cr at 12.3 < log N hi < 12.5 within 
the extrapolated slope f rom ILehner et al. I . 

IPenton et al] (|2004l ') found /3 = 1.65±0.07 (consistent to 
1.5o- with our results) over 12.3 < log Nhi < 14.5 at 0.002 < 
z < 0.069, based on rest equivalent widths and an assumed 
Doppler parameter of fo = 25 km s~^. lDave fc Tripp I (|200ll ) 
used simulations to correct for incompleteness to show no 
break in the H I column density distribution to log Nhi = 
12.6, and calculated a consistent slope to our value. 



4.1.3 Column density distribution 

Obse rvations suggest (e.g. iKim et aL 1 120021: Misa wa et al. I 
l2007f ) and simulations predict fe.g. fPa^ hos ot al. 1 1200^ 
that the column density distribution should steepen at lower 
redshift for log Nhi ^ 13.0— 13.5, which would be consistent 
with a systematic decrease of column densities over time for 
overdensities of a given strength. Our data can help to in- 
vestigate this trend, by uniquely probing the very lowest 
H I column densities in the Lya forest at low z. However, 
given the small sample size and limited number of absorbers 
with log A'^Hi ^ 13.5, we cannot obtain a good log dN/dNni 
slope constraint using just our data, in the framework of the 
commonly adopted parametric form dN/dNni oc Nhi~^ . 
Therefore, we employ additional data from the literature 
to provide a larger comparison sample (larger Aa), with 
the majority of absorbers at higher Nhi to extend the 
column density range. We used the Lya forest sample of 



20 G. M. Williger et al. 



\ 




Figure 5. HI column density distribution, log[/(A'^)] = 
\ogd?N/dXdNHi. We plot the \ogN„i > 13.2, h < 40 km s-\ 
< 40% error in b and Njjj subsample from lLehner et al. 1 1|2007^, 
recalculating the norm alisation to correct an error in their fig. 14 
i Lehner et al. iboOTlJ) . Diamonds: dN/dNui for iLehner et al~1 
I 2007al) Their data are stated complete to log Nhj > 13.2, with 
a turnover at lower Nhi- Large Triangles: 3C 273 data, with the 
same Doppler parameter and H I column density restrictions (11 
absorbers). The vertical bin centre bars for the 12.3 < log Nhi < 
12.6 and 12.8 < log Nhi < 13.2 bins have been slightly off- 
set to the right for clarity. Small Triangle: corrected value for 
12.3 < log Nhi < 12.6 bin (six absorbers), assuming a 50% prob- 
ability of detection, estimated f rom our simulations. Solid line: 
our best fit slope to the complete lLehner et al. I ll2007al ) data only, 
log[/(Af)] = ANhi'^, logA = 12.11+o l°, /3 = 1.793 ± 0.066. 
Dashed lines: la range in slope and normalisation. We find the 
column density distribution at 12.3 < log Nhi < 13.2 consis- 
tent at the ~ 1 — 1.5(T level with an extrapolation from that at 
log Nhi > 13.2. A normalisation error of 0.03 dex is negligible 
and not inclu ded. Dotted line: Column density distribution from 
models of iPaschos et al~l |2Q09), table 4, for 2 = 0.1. A normalisa- 
tion correction of log /(A'^) -1-0.2 (60% increase) would best match 
the data. 



The Low-Redshift Lya Forest toward 3C 273 21 



4.2 Broad lines 

Broad Lya absorbers (BLAs), with 6 > 40 km s~^, have 
been a subject of interest in the past few years (e 



Bowen et al. 2002'; 'Penton et al."2004: Richter et aL^ 

Scmbach ct a l. 2004; Richter et al. 2006a; Rich ter et al. I 



2004 



2006bl : [Wakker fc Savage II2009I . and references therein), be- 



cause they may be the sign of a significant reservoir of 
warm-hot IGM (WHIM) gas. Following the convention of 
iRichter et al. I . we use a sensitivity limit of 



3 X 10^- 
IsjN) 



-cm"^(kms"^)"^ ^ 7.9 x 10^ 



(1) 



using the signal to noise ratio per (their) 10 km s^^ 
resolution element in the above expression. This yields 
log Nh I /b « 10.9 given our minimum signal to noise ra- 
tio per pixel of ~ 21 in the Lya forest, ~ 2 pixel resolution 
elements and 3.2 km s~^ pixel size. Our BLA sample con- 
tains three Lya absorbers in total {z = 0.07374, 0.07429, 
0.08763), using the continuum generated by AUTOVP plus 
manual modifications around emission lines, and the com- 
pound feature detection algorithm. Assigning reliability to 
BLA detections can involve subjectivity, due to the nature 
of echelle spectra and variety of detection and continuum 
fitting techniques. Blending, low S/N ratio, kinematic flows, 
Hubble broadening and continuum uncertainties could be 
responsible for up to ~50% of BLAs. To keep some sub- 
jective standards consistent between different studies, we 
had one co-author (KS), a member of several collabora- 
tions which have studied BLAs, examine the BLA candi- 
dates. He determined the z = 0.087632 absorber to be 
consistent wi th some "reliable" sa mples of recent studies 
of BLAs e.g. iLehner et al. I l|2007al 'l. Using this identifica- 
tion to divide the BLA sample into "total" and "reliable", 
we find dN/dz — 25 ± 15 for the total BLA sample and 
dN/dz = 8.5 ± 8.5 for the reliable BLA, using a pathlength 
of Az = 0.117 (Fig. [6]). These are within la of the values for 
the tot al and high quality samples in fig. 2 of iRichter et al.l 
l|2006ah {dN/dz ~ 35 and 15 for total and reliable systems, 
respectively), which were based on simulated spectra drawn 
from hydrodynamical simulations. 

Given that the 3C 273 sight line is covered by the 
Sloan Digital Sky Survey spectroscop ic galaxy catalogue 
(DR6, lAdelman-McCarthv et al.l [iooi ). we examined the 
galajcy environment of all BLAs and O VI systems beyond 
Virgo (Virgo systems being studied in iTripp et ai~l I2OO2I : 
IStocke et aUl2006l : I Wakker fc Savage 1 12009| . . for example) 
and not in the 3C 273 host cluster. The one secure BLA 
and one O VI system each have several listed galaxies with 
spectroscopic redshifts within 1.5 comoving Mpc impact 
parameter and Av < 200 km s~^. Their H I redshifts of 
z = 0.087632 and 0.090220 correspond to a line of sight 
distance of 9.8 local frame Mpc. 

The z — 0.087632 BLA has 2 galaxies with spectro- 
scopic redshifts from the SDSS within 1.3 local frame Mpc 
perpendicular to the line of sight ia„ plane) and |A ti| < 
75 km s"^ (tt axis) (|Adelman-McCarthv et al.ll2008l ). The 
next nearest galaxy in the plane of the sky, in increas- 
ing velocity space difference, is 2MASX J12271948+0200408 
(cTp = 2.95 local Mpc, z = 0.088579, \Av\ = 261 km s"^). 
The next closest galaxies within 1.6 local Mpc on the sky 
are three a.t z = 0.0902 ± 0.0002, which are close in redshift 
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Figure 6. Cumulative number of BLAs per unit redshift, 
log(iiAf/(i2)BLAi ^ a- function of the absorption s trength , 
log(Ar(cm-2)/6(kms-l)) modelled after IRichter et al.l l l2006al 1. 
The total sample is shown in grey and the high quality sample is 
plott ed in black. Th e diamond shows the value for log(dAf/d2)BLA 
from IRichter et al.l The triangles show the data from this study, 
both the total sample (upper) and the secure sample (lower). Ver- 
tical error bars have been slightly offset for clarity. 



to the O VI system (see below). We chose 10,000 random 
points at 0.085 < z < 0.095 over 1° < (5 < 3° in declina- 
tion and 145° < a < 215° in right ascension, which contains 
1825 galaxies allowing a 1° border. The SDSS spectroscopic 
sample covers 15.0 < r < 17.8 in Petrosian r magnitude, 
which corresponds t o 0.4 — 5L* over 0. 85 < z < 0.095 
(M; = -21 A. -21.6 iBlanton et al. II2OOII : iNakamura et al. I 
I2OO3I ). By using the SDSS as its own control sample, system- 
atic undercounting of dense regions due the maintenance of 
a minimum fiber distance should affect both the observed 
and the control samples similarly. We found that a ran- 
domly placed absorber would have at least 2 SDSS galaxies 
within 1.3 Mpc and |Aw| < 75 km s~^ with a probability of 
P = 0.006 — 0.008, using the above search regions and vary- 
ing the redshift limits by ±0.025. There is a third galaxy 
with a redshift in the NASA Extragalactic Database but 
not in the SDSS, LDTA 139870, which is 2.2 local frame 
Mpc from the line of sight with Av — 7 km s^^; if it were in 
the SDSS, the random probability of finding three galaxies 
within 2.2 Mpc would be the same. 

The z = 0.090220 O VI absorber (logO VI = 13.18 ± 
0.06, iTripp et al.1 12008| ) has 3 SDSS galaxies within 1.6 
local frame Mpc perpendicular to the line of sight and 
|Aii| < 50 km s~^. The next closest galaxies within 2.6 local 
Mpc on the sky are the three a.t z = 0.0876±0.0002 discussed 
with the BLA above, and 2dFGRS N387Z108 (o-p = 2.84 lo- 
cal Mpc, z = 0.0920, |At;| = 489 km s"^). We chose 10,000 
random points as for the z — 0.087632 BLA, and found 
that a randomly placed absorber would have at least 3 
SDSS galaxies within the same region with a probability 
of P = 0.0004 — 0.0015, using the above search regions and 
varying the redshift limits by ±0.025. Properties for galaxies 
near this O VI absorber and the reliable BLA are listed in 
Table ISl 
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Table 5. Galaxies near z 
absorber 



0.087 BLA and z 



0.090 O VI 



name 


position (J2000) 


z 


Av (km s ^) 


Ae' 


Dx 


,.6 


references'^ 








from BLA/O VP 




local Mpc 






SDSS J122906. 84+021348.2 


12:29:06.8 +02:13:48 


0.087361 


-75 


10.7 


1.05 


17.0 


2 


LEDA 139870 


12:28:20.8 +02:21:34 


0.087660 


7 


21.7 


2.14 


18.0 


1,3 


SDSS J122828. 60+021109.3 


12:28:28.6 +02:11:09 


0.087717 


23 


12.4 


1.22 


17.6 


1,3 


2MASX J12285184+0206033 


12:28:51.9 +02:06:03 


0.090046 


-18 


4.7 


0.47 


17.1 


1,2 


2MASX J12280760+0202524 


12:28:07.6 +02:02:52 


0.090282 


47 


14.8 


1.50 


17.3 


4 


SDSS J122806. 93+015959.5 


12:28:06.9 +01:59:59 


0.090333 


61 


15.3 


1.54 


16.7 





" Velocity difference from OVI system is for highest column 

density HI component. 

>> Magnitude from SDSS pR6 llAdeIman-McCa,rthv et al.ll200Sf). 
M* ^ -21.5 at 2 ~ 0.09 l lBlanton et al. Il200ll : iNakamura et al. I 
l2003t) . or r RJ 16.5. 

references: (1) iMorris et al~l lll993l) (2) 



Bahcall fc Bahcall I lll97(t), (31 - 

Tzanavaris. Georgantopoulos fc Georgakakis I l|2006l 'l 



Paturel et aFI 1I2OO0I) . (4) 



The closest three galaxies in velocity space to the z = 
0.087 BLA (including LEDA 139870, in the NASA Extra- 
galactic Database but not in the SDSS spectroscopic cata- 
logue) form a rough rhombus with the 3C 273 line of sight, 
implying that the BLA gas is located in a sheetlike struc- 
ture at least ~ 1.1 x 2 local Mpc in extent. The multi- 
phase O VI absorber at zhi = 0.089898 - 0.09110 is at 
velocity separation Av = 682 km s~^ from the z = 0.087 
BLA. The highest H I column density component of the 
O VI absorber is —18 < Av < 61 km s"'^ from the three 
z ~ 0.09 galaxies, and they trace an elongated forma- 
tion ~ 1.5 local frame Mpc long. The six galaxies are 
moderately luminous {L/L* ~ 0.5 ± 0.2). Based on SDSS 
spectra and colours, and employing SED fitting from PE- 
GASE ijPioc fc Rocca-VolmerangeTll997i ). all six galaxies 
near the BLA and O VI absorber appear to be moderate 
mass (9.8 < logM/Mg < 10.9), old (>10 Gyr) and with 
low star formation rates (0.1 — 1.7Mq yr^^). See Fig. [7]for 
the relative positions of the 3C 273 sight line and galax- 
ies. We also note that there is a QSO/AGN close to that 
redshift, 2MASx J12260241+0046403, which has z = 0.083 
(Av = 1200 km s"^), b j = 18.57 and is 89.3 arcmin or 8.4 
local frame Mpc away I Hao et al. II2005I ). which may be in 
a filament with the BLA. The fact that the two absorbers 
are separated by ~ 700 km s~^ and that each is on the scale 
of 1-2 Mpc from moderately bright galaxies may suggest 
a correlated physical environment for both absorber-galaxy 
associations. 

The BLA and O Vl-galaxy distances of 1-2 Mpc 
are consistent with the correlation lengths for such ab- 
sorbers with ^ 0.1 — l.OL * gala xies found from sim- 
ulations by iGangulv et al. I (|2008l ) but are larger than 
the 1 00-300 kpc values found by lOppenheimer fc Dave I 
l|2009l ). Observationally, the velocity differences between 
the galaxies and BLA+0 VI absorber are well within the 
~ ±300 — 500 km s~^ differ ences used for much of the 
IWakker fc Savage I (|2009l ) and IStocke et al~l l|2006l l analy- 
ses, almost to the point that any putative absorber-galaxy 
structures for our two examples appear remarkably coher- 
ent over Mpc scales. The galaxy-absorber distances are a 
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Figure 7. Schematic of 3C 273 (cross) with galaxies within 
\Av\ < 200 km s-^ of the z = 0.087632 reliable BLA (diamonds) 
and \Av\ < 50 km s'^ of the z = 0.090220 O VI absorber (tri- 
angles). At 2 = 0.09 in our adopted cosmology, 0.1° on the sky 
corresponds to 0.60 local-frame Mpc. 



factor of several larger than the few x 100 kpc distances 
found in a l arge, z < 0.1 5 gala xy- O VI system correla- 
tion study bv lStocke et al. I ([2006), and also are large com- 
pared to the 350-450 kpc range for O VI absorbe r s with 
L > O.I L* galaxies found by IWakker fc Savage I (jioogh 
at z ~ 0. IWakker fc Savage I iProchaska et al. I (|200(: ) and 
ICooksev et al. ( 20081 ) found significant variations in the 
galaxy environment of O VI absorbers. The proximity of the 
secure BLA and O VI absorber with each other and with lu- 
minous galaxies is consistent with the observation of analo- 
gous z ~ 0.06 — 0.08 O Vl-galaxy correlations in large- scale 
filaments (Tripp ct al. 2006) and the WHIM ori gin of B LA 
systems (e.g. iRichter et al.ll2006al : iRichter et al. ir2006d ). 

Finally, we compared the general correlation between 
local galaxy density and in tegrated Nhi column density for 
the 3C 273 sight line as in lBowen et al. 

] (|2002l ). and made 

a Spearman rank correlation and Kendall r correlation test 
similar to the one in Paper I. We used 43 galaxies from the 
NASA Extragalactic Database within 2° and 2 local frame 
Mpc of the 3C 273 sight line at 0.020 < z < 0.139, where our 
Lya detection probability is > 80%. The galaxy sample is 
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inhomogeneous and dominated by SDSS entries, but with no 
more significant selection effects tlian the sample used in Pa- 
per I. We do not probe further out in angular distance due to 
incomplete coverage in the SDSS DR6. We find no evidence 
of a correlation between galaxies and column density using 
absorbers with log Nhi > 12.5. The Spearman rank coeffi- 
cient for redshift bin Az = 0.003 {Av ~ 800-900 km s~^) is 
p = 0.283 with probability pp — 0.085, and the correspond- 
ing Kendall rank correlation is r = 0.250 with pr — 0.027. 
Results are insensitive to binning on larger scales, and to 
lowering the redshift cutoff to z = 0.11, because there is only 
one additional galaxy at z — 0.138 below the end of the Lya 
forest sensitivity limit. Smaller scale bins are dominated by 
zero counts of galaxies and H I column densities. This lack 
of a correlation contrasts with the weak correlation between 
galaxy number density and summed Nhi column density on 
4000-6000 km s~^ scales detected toward PKS 0405-123 in 
Paper I. 

The 3C 273 sight fine should be probed to <0.1L* 
up to the QSO redshift to determine whether less lumi- 
nous galaxies at smaller distances can be associated with the 
O VI and BLA absorption. Larger cross-correlation samples 
of BLAs+galaxies and BLAs-|-0 VI systems are needed to 
quantify further any such trends. 



4.3 Clustering and voids 

4-3.1 Clustering 

Clustering in the low z Lya forest has been studied in a 
number of cases (e.g. Ijanknecht. Baade Sz Reimers 1 12002| . 
and references therein), and weak clustering has been indi- 
cated on velocity scales of At; < 500 km s^^. We found a 
clustering signal on the Av < 250 km s~^ scale with STIS 
E140M data covering z < 0.4 in a sample of 60 Lya absorbers 
with logNHi > 13.3 toward PKS 0405-123 (Paper I). The 
clustering toward 3C 273 is expected to be weaker than for 
PKS 0405-123 due to the lower Nhi threshold. In a manner 
similar to Paper I, we established a minimum line separation 
of Av = 31 km s~^, based on the observed line separation 
distribution in the actual data, and made 10000 random 
simulations of the 3C 273 line list for c omparison to the 
data, assuming dN/dz oc (I + 2:)'', 7 = 0.26 lWevmann et al. I 
l|l998ll . 



We find clustering at the 3.0a significance level in the 
two point correlation function ^(Av) for the 21 main sam- 
ple absorbers with log Nhi > 12.5 at 0.020 < z < 0.139, 
on the scale of Av < 1000 km s~^. The signal increases 
in significance to 3.4 a, corresponding to a strength of 
^{Av < 1000) = 1.0 for log Nhi > 12.6. Poisson counts 
of the pairs in the bin, which are only an approximation 
as the pairs between bins are correlated, would indicate 
la errors of ~ 40% . The simulations give a 99% con- 
fidence upper limit of ^{Av < 1000) = 1.1. Based on 
the 10000 simulated line lists, the probability of the sig- 
nal being matched or exceeded is P = 2.8%. The signal 
then diminishes in strength at even higher column den- 
sity thresholds due to the declining sample size. Using the 
same hydrodynamical simulation s for comparison as Paper I 
(|Dave. Katz. fc Weinberg1l2003l '). we find the measured clus- 
tering limits consistent wit h the predictions (Fig. [8]). The 
r-step correlation function (|Kim et al. Il2003 ) shows no sig- 



> 



c3 



?3 




2 4 6 

AV (/lOOO km s"^) 



Figure 8. Two point correlation function for the Lyo sample 
(log Nhi ^ 12.6) with Av = 1000 km binning, which shows 
the strongest correlation in our sample, (3.4 a). Shaded regions 
denote 68, 95, 99% confidence limits from 10* Monte Carlo sim- 
ulations. Simulated line lists have been filtered to eliminate pairs 
with velocity splittings of Av < 31 km s~^, as in the data. Error 
bars show example la Poissonian erro rs. The triangles show ex- 
pected values from the simulations of lDave. Katz, &: Weinberg I 
||2003|) for the same column density threshold; the simulation box 
size artificially lowers the value of § in the higher two velocity 
bins. 



nificant signal. The function appears sensitive to rare high 
column density systems in an otherwise small sample dom- 
inated by weak lines, as the only feature of note was con- 
tamination at the All ~ 70 — 90 km s"'^ scale for r > 0.4 
from the z = 0.066548 absorber with log Nhi = 14.08. The 
restricted Lya sample with 6 < 40 km s'^ and errors in 
b and Nhi of < 40% only contains 16 systems, and is too 
small for useful clustering analysis. 

Our weak signal in the low z log Nhi ^ 12.5 — 12.6 Lya 
forest two point correlation function for Ai; < 1000 km s~^ 
is consistent with the weak clustering predicted by numerical 
models. However, it is on a larger scale than the signal found 
toward PKS 0405-123 {Av < 250 km s~ ^) ior log Nhi ^ 13.3 
On a similar Aw = 1000 km s~^ scale, iBowen et al. I (I2OO2I ) 
found that Nhi column densities summed on 1000 km s~^ 
scales correlate well with the volume density of Mb < —17.5 
galaxies, which was reminiscent of the effect we found in 
Paper I, albeit on 4000-6000 km s~^ scales. 

There have been no comparable clustering studies in 
the low redshift Lya forest to column density thresholds as 
low as log Nhi ~ 12.5. Weak clustering has been found in 
the low z Lya forest for stronger abs orbers on scales of up 
to 500 km s~ MKirkman et aL [(|2007l l noted non-zero auto- 
correlations in 74 FOS spectra with 230 km s"'^ resolution 
at z < 1.6 out to Av < 500 km s~^ using pixel opaci- 
ties, which were strongly peaked at Av < 100 km s~^ for 
z < 0.5 and more broadly peaked at 0.5 < z < 1.5. Their 
column density limits arise from sight lines with a variety 
of sensitivities, but are pres umedly ^,0.5 — 1 . dex higher 
in column density than here. Ijanknecht et al."l (|2006l ) noted 
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weak (^(Aw) ~ 0.1 — 0.4), marginal signals on scales ol 100 < 
An < 200 km s ~^ and 1000 < Av < 2000 km s~^ for 12.9 < 
log Nhi < 14.0. |Penton et al.l (|2004l l determined a two point 
correlation function signal of ^(Ad < 190 kms~^) ~ 3.3 at 
the 4.5(T significance level and ^(Ai; < 260kms~^) ~ 2.8 
at 5.6(T significance for rest equivalent width Wq > 65 mA 
{log Nhi >IS.1 for 6 = 25 km s"^). They also found a 3cr 
excess at 260 < Av < 680 km s~^. However, they found 
no significant correlation for weak absorbers (Wo < 65 mA) 
with either themselves or with stronger absorbers, which 
is consistent with our results. Our 99% confidence limit of 
5(Aii) = 1.1 at Av < 1000 km s~^ is at least crudely com- 
parable if not a slightly lower constraint than their four la 
upper limits of ^(Aw)<2 over a similar velocity range. A 
sample size an order of magnitude larger (containing 100 x 
more pairs and roughly reducing bin errors by a factor of 3) 
would make either more interesting clustering constraints 
for weak systems or detect a correlation strength predicted 
in our comparison models. COS data should address this 
question. 



4.3.2 Voids 

As in Paper I, we searched for regions in velocity space 
empty of Lya systems above a series of H I column den- 
sity thresholds, and compared them with the frequency 
of similarly-sized gaps in randomised line lists having the 
same statistical properties as the observations. The most 
significant gap we find has a line of sight extent of An = 
8100 km s"^ (120 comoving Mpc) at 0.0901 < z < 0.1200 
for log Nhi > 12.6. The probability of such a gap being 
matched or exceed in a sample with the same number of 
Lya lines is P = 0.006, using simulated line lists having a 
Poisson mean of 21. The void persists at more marginal sig- 
nificance for log Nhi > 12.7, with P — 0.022, and for higher 
Nhi thresholds the sample becomes too small to be statis- 
tically useful. For log Nhi > 12.5, the most significant void 
is over 0.0072 < z < 0.0262, with P = 0.023. 

We can use the SPSS D R6 spectroscopic sample 
lAdelman-McCarthv et al] |2008| ) as a rough indicator of 
galaxy overdensity, although by 2 ~ 0.12 it is only prob- 
ing bright galaxies. We count 0,3,9,14 galaxies in the SDSS 
DR6 within 1,2,3,4 local-frame Mpc respectively in the line 
of sight at 0.0901 < z < 0.1200. We made counts of galaxies 
within a grid of 1-4 local frame Mpc around 10000 random 
locations in the same redshift interval, over 140° < a < 225° 
and —0.25 < 5 < 3.0°, and found a minimum random 
probability (for galaxies out to a maximum radius of 
1.48 Mpc) of P = 0.88 for the Lya forest gap to be in 
a region up to similarly as dense. The max;imum proba- 
bility is P = 0.988 - 0.989 for radii r = 3, 4 local Mpc. 
The galaxy environment therefore does not appear signifi- 
cantly underdense, and could in fact be overdense. A com- 
parison with the interval 0.08 < z < 0.09, which is com- 
paratively rich in that it contains 3 absorbers and 0,2,4,6 
galaxies within r < 1,2,3,4 local frame Mpc, respectively, 
reveals no significant over or underdensity. (P ranges from 
0.771-0.922.) Deeper galaxy surveys may reveal more about 
the relationship between the weakest absorbers and nearby 
galaxies. We do note from the NASA Extragalactic Database 
(NED) that a QSO ~ 1° from 3C 273, RX J1230.8+0115 



has V=14.4, rest-frame Lyman limit flux F ° ~ 5 ± 1 x 

■' 1018A 

10"^* erg cm^ s"^ A"^ (ProS program 1099001, PI Sem- 
bach), z = 0.117 and is 54. 3 arcmin or 6.9 local frame Mp c 
from the 3C 273 sight line l|Read. Miller fc Hasing"er1ll998l '). 
For comparison, the QSO SDSS J030435.32-000251.0, which 
was suggested as causi ng the z = 3.05 H e II opacity gap to- 
ward Q0302-003 (Jakobsen et al. 112003 : Heap ct al. 2000), 
is only 3.0 local frame Mpc from the line of site and is likely 
brighter by several tenths of a dex than RX J1230. 8-1-0115 
(u — 22.35; a full comparison is beyond the scope of this 
work). The z = 0.1200 O VI absorber toward 3C 273 is 
only Av ~ 800 km s"^ from the RX J1230.8+0115 redshift, 
and both objects may be located in the same filamentary 
structure. 

Voids in the Lya forest are rare, and most of the work 
in the field has been done at high redshift. Few Lya forest 
void studies have been done at z < 1.6, but they are valu- 
able as comparisons to the galaxy- absorber relation with 
respect to the inverse questio n of Lya absorbers in galaxy 
voids fe.g. IStocke et aLll2007^ . In Paper I, we found toward 
PKS 0405-123 gaps at 0.0320 < z < 0.0814 of 206 comoving 
Mpc as defined by log Nhi > 13.3 absorbers with probabil- 
ity of random occurrence in shuffled data of P = 0.0004, at 
0.1030 < z < 0.1310 for log Nhi > 13.2 (113 comoving Mpc, 
P = 0.007) and at 0.0320 < z < 0.0590 for log Nhi > 13.1 
(also 113 comoving Mpc, P = 0.003). There are several 
higher redshift voids {z ~ 1.6 - 3.2) on the 40 - 60 Mpc 
scale for log Nhi ~ 13.4 — 13.5, which we discussed in Pa- 
per I. Given the expected column density vs. density per- 
turbation evolution between z ~ 3 and z ~ 0.2, the column 
density levels for which we find our low redshift void(s) to- 
ward 3C 273 (12.5 < log Nhi < 12.7) are comparable to 
log Nhi ~ 14.0 at jz = 3 (§ refsubsubsec:redshiftdensity). 
Our 3C 273 data are thus ~ 0.5 dex too shallow to probe 
similar density perturbations as used to define the gaps at 
z ~ 2 — 3. Therefore, we do not consider it unusual to find 
a slightly larger candidate gap than reported at higher red- 
shift, given the difference in column density sensitivity. We 
do not find the 0.0901 < z < 0.1200 region within radius 
r < 4 local frame Mpc of the 3C 273 sight line to be of 
unusually low galaxy density, based on the SDSS DR6 spec- 
troscopic survey. A more definitive measure of the galaxy 
environment would require a more sensitive survey. 

Although the 0.0901 < z < 0.1200 gap toward 3C 273 
is coincident with RX J1230.8-I-0115 {AO = 54.3 arcmin 
or 6.9 local frame Mpc), recent analyses of any trans- 
verse proximity effect on smaller scales of ~ 1 — 5 Mpc 
have not indicated significant thinning of the Lya for- 
est which could produce voids. Possible explanations for 
the lack of any detections include anistropic emission (e.g. 
Hcnnawi & Procliaska 2007), or the masking of any trans- 
verse proximity e ffect by enhanced density env ironments 
around QSOs (,Goncalves. Steidel fc Pettini II2OO8I '). The cor- 
relation length between QSOs and optically thick H I ab- 
sorbers off t he line of sight is 9.2'^ ]'^h~^ comoving Mpc 
at z ~ 2.5 (|Hennawi fc Prochaska \ ). which is three times 
stronger than the QSO-Lyman break galaxy correlation, 
would be on the order of the observed distance between 
RX J1230.8+0115 and the z = 0.1200 O VI absorber to- 
ward 3C 273. Both objects could therefore mark the same 
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structure. A deep galaxy survey toward the region would 

conf irm this. 

IShang. Crotts fc Hairiianl (j2007l ) searched the SDSS for 
gaps using a pixel opacity technique with limits of transmit- 
ted flux F = 0.6 - 0.8 at 2.8 < 2 < 3.9, which probes rest 
equivalent widths corresponding to log Nhi ~ 13.0 — 13.5 
for 18 < & < 45 km s~^, or overdensities of logp/p — —0.2. 
They found void sizes of ~ 3 — 35 comoving Mpc, which 
were in agreement with the standard ACDM model for Lya 
forest formation. As our Lya absorber sam ple does not 
probe down to the same overdensity level as the lShang et al.l 
search, our results are not inconsistent with theirs. Given 
the probability of random occurrence of the largest gap we 
identified, the current paucity of low redshift data sensitive 
to log Nhi ~ 12.5 a nd the lack of data which probe the 
same overdensities as IShang et al.l . future studies with COS 
should be used to confirm the frequency of large voids. 



4.4 The Virgo cluster absorber at z — 0.00530 

iTripp et al. I (j2002) established metallicity limits and corre- 
sponding physical conditions for the two Virgo cluster ab- 
sorbers, including the one at z = 0.00530. We took the op- 
portunity to re-fit the components. We used nine Lyman se- 
ries lines through LyK (919A), and included data from STIS 
and the revised reductions from the LiFlA, LiF2B, SiC2A 
and SiClB F(75'£' channel segments. (The SiC segments cov- 
ering Ly/3 were of such low signal that they did not help to 
give useful constraints.) 

We took advantage of the clear Si III feature to set 
the redshift for the main H I component. Two components 
provide a satisfactory fit to the data (Fig. [9]). The reduced 
~ 1.19, though the formal probability of a fit to the 
data is still < 1%. One component makes a significantly 
poorer fit, and three are not substantially better (reduced 

~ 1.234). A four component fit is marginally worse 
than the two component fit (reduced = 1.209), but suf- 
fers from increasing uncertainties in the H I column den- 
sities. We therefore adopt a two component fit with a to- 
tal H I column density log Nhi = 15.70 ± 0.04 {la error). 
This is within the la error of the one component fit from 
the curve of growth measurement of log Nhi = 15.77to' ;[o 
excluding Lya, which has the mos t profound blen d ing e f- 
fect by the second component, bv ISembach et alTI l|200 j ). 
iTripp et al.l (12002) ado pted log Nhi = 15.85lg jg, which 
was the Sembach et al. I figure (including Lya in the curve 
of growth fit) f or their metallicity studies. The conclusions 
of lTripp et al. I about the nature of the z — 0.0053 absorber 
should not be significantly affected. It is possible that the 
component structure is more complex, but higher signal and 
resolution far-UV data would be necessary to make a signifi- 
cant improvement, which must wait for a successor to FUSE. 
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Figure 9. Two component fit to ^ = 0.00530 Virgo absorber, 
including Si III 1206 and C II 1334. Green - total profile. Red - 
higher H I component, at the metal z. Blue - blue component, 
H I only. Lilac - blends: solid - H2; dashed — unknown (UNK). 
Note that the colour sequence from top to underneath is lilac- 
blue-red-green where there is more than one component (i.e. H I, 
not metals) so sometimes green especially is obliterated by a com- 
ponent. 
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5 CONCLUSIONS 

We have performed an analysis of the Lya forest toward 
3C 273 based on STIS E140M and a re-processing and re- 
duction of FUSE spectra, and find the foUowing. 

1. We present STIS EI40M echelle data for 3C 273, and 
performed profile fits or measured apparent optical depths 
for all of the detected Lya forest over Q < z < 0.158, plus 
intervening metal and Galactic metal absorption systems. 
We analysed simulated spectra to determine our sensitivity 
to line detections and resolution of close pairs, and probe 
to log Nhi ~ 12.5 over most of the spectral range. Our 
main sample consists of 21 absorbers to log Nhi ~ 12.5 over 
0.020 < 2 < 0.139. 

2. The redshift density for the absorbers with column 
densities log Nhi > 14.0 is consistent with previous, lower 
resolution studies, though based on only 4 absorbers. For ab- 
sorbers with 13.1 < log Nhi < 14.0, the redshift density is 
~ 1.5(T lower than the mean of comparable low redshift stud- 
ies at the same resolution. Absorbers with log Nhi > 12.5 
and a subset using the s ame error-limi t ed, b < 40 km s~^ 
sample criteria as used bv lLehner et al. 1 (|2007al ) both have a 
simila r redshift density to p redic tions from n umerical mod- 
els bv lPaschos et al. I l|2009l ) and lDave etHI (fl999. ). 

3. The Doppler parameter distribution has a mean, me- 
dian and standard deviation of 28, 27, 13 km s~^, respec- 
tively. There is no significant evidence that line blending 
affects the distribution. If we choose a restricted sample 
with errors < 40% in b and log Nhi and remove absorbers 
with fe > 40 km s~^, the mean and standard deviation 
are 24, 24 and 8 km s~^, respectively, which is consistent 
both w ith the main Lya system sample from Leh ner et al. I 
l|2007al ) and also with their similarly restricted data set for 

6 < 40 km s"\ 

4. The H I column density distribution to log Nhi = 
12.5 is consistent with a power law fit to a s ample 7 times 
larger covering 13.2 < logA''£f/ < 16.5 from iLehner et al. I 
l|2007al ). If we correct our sample for the probability of de- 
tection at 12.3 < log Nhi < 12.5, there is no evidence for a 
significant break in the power law fit to log Nhi = 12.3. 

5. We find a total of 3 broad Lya absorbers (BLAs) 
with 6 > 40 km s~^ and log Nhi /b ~ 10.9, one of which is 
deemed reliable. The redshift densities are c onsistent with 
results from hydrodynamical simulations by iRichter et al.l 
l|2006al 'l. The reliable BLA at z = 0.087632 is « 700 km s"^ 
from an O VI absorber in velocity space. The z = 0.0902 
O VI absorber is in a high galaxy density environment with 
a probability of being in an equal or richer environment of 
P « 0.001, based on galaxy counts in the SDSS DR6. The 
reliable BLA is in a marginally overdense environment as 
well {P ~ 0.007), and both may be embedded in the same 
structure. 

6. We find evidence for weak clustering in the Lya for- 
est at 3.4cr significance for logA''(H I) > 12.6 from the two 
point correlation function on a scale of Aw < 1000 km s~^. 
The limits obtained are consistent with the clustering ex- 
pected at that column density limit from hydrodynamical 
simulations. 

7. We find evidence for a void in the Lya forest for 
log Nhi > 12.6 at 0.0901 < z < 0.1200, with a random prob- 
ability of occurrence P — 0.006. The redshift interval does 
not cover an unusually low galaxy density environment. The 



relation between weak absorbers and galaxy environment as- 
sociated with the 3C 273 and other low z sight lines should 
be addressed with deeper galaxy surveys than the SDSS. 

8. A re-analysis of the z — 0.0053 Virgo Cluster metal 
absorption system using both FUSE and STIS data reveals 
velocity structure in the H I absorption which is best fitted 
by 2 components. The total H I column density is within 2a 
of that used for a previous analysis of the physical properties 
of the system, whose results are not significantly changed by 
the adoption of a relatively weak second absorber. 

3c 273 provides one of the best probes to very low col- 
umn density limits at low redshift, and thus a unique probe 
of halo gas around galaxies in the local universe. Further 
studies with COS and a deep, complete galaxy survey will 
no doubt reveal further details of the galaxy-absorber rela- 
tion. 
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APPENDIX A: A COMPARISON WITH THE 
LINE LIST FROM DANFORTH & SHULL (2008) 

iDanforth fc ShulTI (|2008l ) published a study of a large sam- 
ple of QSO absorbers, including the line list for a data set 
analysed in detail in this work. There a re a number of dif- 
ferences in the reduction and analysis. IDanforth fc ShuUI 
used Calfuse 2.4 to reduce the FUSE data, whereas we 
used Calfuse 3.1.3 (|Dixon et al. 1 120071) . which has many 
changes, including a significant im provement in the w ave- 
length calibration from version 2.4. IDanforth fc ShulTI nor- 
malised the data in 10 A segments centred on IGM absorbers 
and binned the FUSE and smoothed the STIS data by 3 pix- 
els each, whereas we used a mainly automated continuum- 
fitter with manual adj ustments around em ission lines and 
did not bin or smooth. IDanforth fc Shull | used a combina - 
tion of apparent optical depth ( Sembach fc Savage II 199*3 ) . 
curve of growth concordance plots and Voigt profile fitting, 
whereas we mainly used Voigt profile fits and only used ap- 
parent optical de pth to establish low er limits for saturated 
metal transitions. IDanforth fc ShulTI tended to treat blends 
which they could not separate unambiguously as single ab- 
sorbers, whereas our normal procedure is to add components 
when the existing Voigt profile fit has a probability of match- 
ing the data of P < 0.01. 

We matched absorbers from the IDanforth fc Shull I list 
with ours, and combined apparently blended absorbers in 
our line list in velocity space (up to Aw ~ 70 km s~^) by 
summing their H I column densities. We display the com- 
parison for 20 systems in common in Fig. lAll For the plot we 
use the errors in log Nhi for our fits, and also the Doppler 
value and corresponding errors, for the highest Nhi compo- 
nent. The H I column densities show a good match within 
la errors. The Doppler parameters show more scatter, which 
likely arises from the finer structure decomposition in this 
study, but are still largely consistent within la errors. 



Absorbers in one list but not the other, or potentially 
suspect absorbers, are as follows. 

z = 0.0026 30 : This is a broad line which 
IDanforth fc ShulTI do not list. We need this component for a 
proper fit to the Galactic Lya profile, which we fitted simul- 
taneously with the IGM components in the wing. It is below 
the minimum redshift for any of our statistical analyses. 

z = 0.007 588 : This feature is included by 
IDanforth fc ShulTI as well. It is in the wing of Galactic H I, 
and shows a double minimum. The complex has ~ 3(t sig- 
nificance overall. We include it due to its contribution to a 
fit to the Galactic H I complex, but note that the profile is 
atypical and may indicate a blend of two weaker features. 

z = 0.054 399 : This is a weak line which 
IDanforth fc ShulTI do not list. It is below our 80% proba- 
bility of detection threshold of log Nhi ~ 12.5, but is used 
for the H I column density extrapolation to log Nhi = 12.3. 

z = 0.060 51 : This is a weak line from 
IDanforth fc Shull I which we do not list due to its low signif- 
icance. It contains a small noise spike halfway through the 
feature, and has a square profile and flat minimum. 

2 — 0.06 4078 : This is a s hallow, moderately broad 
feature which lOanforth fc Shull 1 do not list. We find satis- 
factory constraints on the column density and Doppler pa- 
rameter such that we use it in our statistical samples. 

z = 0.07373 8 and 0.074286 : T hese are broad, shal- 
low feature which lOanforth fc Shull I do not list. Their high 
Doppler parameters of 6 = 64 ± 15 km s~^ and b — 
51 ± 14 km s~^ make us omit them fr om our restricted sam - 
ple analysis, to be consistent with the lLehner et al."1 (|2007al ) 
6 < 40 km s~^ sample criteria. The z = 0.073738 system is 
unlikely to be Galactic O I*, since its redshifted corrected 
position is just outside the O I 1032 profile. There may be a 
better case for OI** being the z = 0.07429 feature, since it 
falls in the ground state profile, and there is possibly O I* 
in the red wing of Si III 1304. However, neither agree par- 
ticularly well in position with C II*, so unless they would be 
at significantly different velocities (which is not suggested 
by the relative positions of the C II/ O I centroids) then we 
doubt that the O I*, 01** identifications apply. 

z = 085882 : This is a broad, shallow feature which 
IDanforth fc ShulTI do not list. It is below our 80% probabil- 
ity of detection threshold of log Nhi ~ 12.5, but is used for 
the H I column density extrapolation to log Nhi = 12.3. 

2 — 0.1 09141 : This is a shallow, round-bottomed 
feature which IDanforth fc ShulTI do not list. We use it in 
our extrapolation of the H I column density distribution to 
log Nhi ~ 12.3, but in no other statistical sample. 

2 = 0.13943 Danforth et al. list it, but we find it less 
plausible. It is assymetrical and contains a noise spike. The 
significance is ~ 3o", and depends strongly on the integration 
bounds. 

2 = 0.141676 and z = 0.150204 : These are broad, shal- 
low features which lOanforth fc Shull 1 do not list. They have 
~ 3ct and 4.5cr significance, respectively. They fall on Lya 
emission and within our definition of the proximity effect 
zone {Av < 5000 km s~^), and are also below our 80% prob- 
ability detection limit. We do not use them for any analysis. 

2 — 0.157787 : This is a dist inct absorber with a hint 
of an extended red wing, noted in lTripp et al. I (|2008l ) in as- 
socia tion with O VI. We fitted the feature with one compo- 
nent. iDanfort^^^^ShuU] do not list it, presumably because 
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it falls within their definition of the proximity effect zone 
(Av < 1500 km s"^). 
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Figure Al. Comparison of H I column densi ty and Doppler pa- 
ramet er for 20 absorbers in common between IPanfprth fc ShuU I 
l|2008ll and this work, witli la error bars. Column densities have 
been summed in the case of fitted complexes in this work span- 
ning Ai; < 70 km s~^. However, column density errors for this 
study have been retained for the highest Nni component, as are 
the Doppler parameters and errors. Fits are generally consistent 
with in la errors. Differ ences are qualitatively understandable in 
that iDanforth fc ShuU I have slightly higher column densities and 
Doppler parameters due to their not using as fine component 
structure as in this work. 
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